Global Diffusion in a Realistic Three-Dimensional Time-Dependent Nonturbulent Fluid Flow 
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We introduce and study the first model of an experimentally realizable three-dimensional time-dependent 
nonturbulent fluid flow to display the phenomenon of global diffusion of passive-scalar particles at arbitrarily 
small values of the nonintegrable perturbation. This type of chaotic advection, termed resonance-induced 
diffusion, is generic for a large class of flows. 
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Dynamical systems that arise in problems of particle diffu- 
sion [1] in incompressible and nonturbulent fluid flows, apart 
from being of theoretical interest, hold much relevance for 
technological applications. Properties of emulsions, disper- 
sion of contaminants in the atmosphere and ocean, sedimen- 
tation, and mixing, are just a few examples. In compari- 
son with the present state of knowledge of chaotic advec- 
tion, mixing and transport in two-dimensional [2] or three- 
dimensional time-independent [3] flows, little is known for 
three-dimensional time-dependent flows. 

The strobed dynamics of fluid parcels — so called passive 
scalars — in three-dimensional time-dependent incompress- 
ible fluid flow is qualitatively equivalent to the iteration of a 
three-dimensional volume-preserving map. A few years ago 
we began to investigate nearly-integrable classes of such maps 
— which we have termed Liouvillian maps — in search of 
possible generic features [4]. We found that we could char- 
acterize Liouvillian maps in a similar manner to Hamiltonian 
dynamical systems, by using the number of slow action and 
fast angle variables that, in the integrable limit, remain in- 
variant and rotate uniformly, respectively. The cases having 
one or two of such actions were found to be particularly in- 
teresting. For the one-action case, which has been shown to 
be generic for a class of flows [5], a KAM-type theorem ex- 
ists [6], and with it the associated barriers to global transport 
[4]. On the other hand, the case of two actions, which is also 
generic for a class of flows, displays in Liouvillian maps a 
new phenomenon of resonance-induced diffusion leading to 
global transport throughout phase space [4]. However, de- 
spite this theoretical progress, until now no realistic flow — 
in the sense of being at least an approximate solution of the 
Navier-Stokes equations for a realizable experiment — has 
been observed to show the above properties. The purpose of 
this Letter is to fill the gap for the two-action case by intro- 
ducing the first engineerable flow demonstrating the presence 
of resonance-induced diffusion in real fluid flows. 

We begin by briefly illustrating the phenomenon for a map, 
introducing at the same time the adiabatic invariants and res- 
onant surfaces we encounter in two-action flows. The map 

I[ = h+ sFrih, 6), I' 2 = h + sF 2 (F l ,6), 

0' = + u(I' l ,I' 2 ), (1) 



represents a small perturbation of an integrable Liouvillian 
map having two action and one angle variables [7]: for van- 
ishing e, the two action variables I\ and I 2 remain invariant 
while the angle 6 rotates with a constant angular frequency 
that depends only on the actions. When e is nonzero but small 
the action variables drift slowly compared to the angle. This 
separation of scales allows us to average over 6 leading to an 
adiabatic description of the motion. Before the actions change 
appreciably, the angle is able to traverse a large sample of its 
domain. Hence the evolution of the actions is sensitive only 
to the average value of the angle: the action equations can be 
averaged over the angle, and so become decoupled from the 
angle equation to yield an area-preserving map 

I[=h+ eFi (I 2 ) , I' 2 = h + eF 2 (,I[ ) , (2) 

where the bar represents the average. Being a small pertur- 
bation of the two-dimensional identity, the map of the action 
plane can be well approximated by a two-dimensional au- 
tonomous Hamiltonian flow 

FlH - FlH 

h = Hh) = Wi , h = F 2 (U) = - — . (3) 

One is then led to the conclusion that the dynamics of the map 
in Eq. (1) occurs on invariant surfaces that are the product of 
the almost uniform angular motion of 9, and the level curves 
of the Hamiltonian 

H(h ,h)= [ 2 F\ (F)dl - [ ' F 2 (F)dI = (]. (4) 
■Jo Jo 

However, this analysis breaks down when the angular vari- 
able does not sample its domain of variation uniformly. This 
failure is bound to occur when the rotation is resonant; when 
nuj{I\, I 2 ) = 2irk, where k and n are integers. The family 
of curves in the action plane defined by Eq. (4) is generically 
transversal to the curves on which these resonances occur. 
Although the resonances are dense in the action plane, they 
have a hierarchical structure in which the relative strength of 
different orders k/n of resonances is governed by the Fourier 
expansion of Eq. (1), such that in general the lowest-order res- 
onances have the smallest denominators n. Hence one might 
expect the adiabatic approximation to provide a good descrip- 
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tion of the action motion, except at the intersections with the 
lowest-order resonances. This is depicted in Fig. l(a)forapar- 
ticular two-action Liouvillian map. Close to resonances the 
trajectory oscillates wildly between different invariant curves 
H = j3. In Fig. 1(b) we have plotted the time evolution of 
the adiabatic invariant H for this map to show that it remains 
nearly constant almost all the time, but jumps chaotically from 
one invariant curve to another at each intersection. 
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FIG. 1 . (a) Projection of a trajectory in a two-action Liouvillian 
map onto the plane z = 0. Themapisx' = £+0.001 (sin 2+2 cos y), 
y' = y + 0.001 (1.5 sin x + 2.5 cos z), z = z + 4(cosj/' + sinx'). 
The diagonal dashed line shows the location of the lowest-order 0/ 1 
resonance, (b) The time evolution of the corresponding adiabatic 
invariant H = 2 sin y + 1.5 cos x . 

The study of two-action properties in a real fluid flow is 
made difficult by the fact that very few three-dimensional 
time-dependent fluid flows are analytically tractable. In the 
following we introduce a fluid flow that can be solved an- 
alytically as a perturbation expansion of the Navier-Stokes 
equations and that exhibits the phenomenon just illustrated. 



The flow we have chosen to consider is incompressible flow at 
low Reynolds numbers between two concentric spheres which 
rotate with different angular velocities about a common axis 
that switches in turn between two different directions separated 
by an angle a. In the absence of time dependence introduced 
by the axis switching, the flow consists of a primary spherical 
Couette flow about the rotation axis, superposed with an or- 
thogonal secondary flow in the meridian plane (r, 0) [8]. Ow- 
ing to axial symmetry, none of the flow components depend 
on the azimuthal coordinate <f>. The secondary flow, which is 
made up of one or two Taylor vortices in each hemisphere, 
depending on the flow parameters [9], is then two dimensional 
and described by a stream function t/> that can be computed 
perturb atively in the Reynolds number Re. The velocity field 
up to first order in the Reynolds number reads 

1 dip 1 dip ■ 

Vr = 2 • a ~*E = r ' V0 = ^~E~s~ = r9 > 

r l sin o0 r sin or 

= (a\r + ^jj sin 6*= r sin 6* </>, (5) 
where ip is given by 

ip = Re + A 2 + A 3 r 3 + A 4 r 5 

+ ^ -a,r 2 )) sin 2 cos 0, (6) 

and ai, a 2 , A\, A 2 , A 3 , A 4 , are constants dependent on the 
ratios of the radii and the angular velocities of the spheres [9]. 
Notice that for Re <C 1 the secondary flow (v r , vg) is very 
much slower than the primary flow v$ . 

The perturbation of this flow by periodic axis switching 
introduces time dependence at the same time as true three 
dimensionality by coupling together the velocity components 
of Eq. (5). Since Re < 1, we can assume that fluid inertia 
is not important, such that trajectories of passive scalars will 
be a piecewise juxtaposition of the steady flow about each 
axis. Notice that the r and coordinates are always coupled 
by the secondary flow except when the Reynolds number is 
precisely zero. Furthermore, adding a second rotation axis 
leaves the equation describing the evolution of r unchanged, 
but introduces a coupling between and <f> which is of 0(a) for 
small axis separation a. All this implies that, although we do 
not have direct access to the three-dimensional integrable case, 
we can however intuit that it should exist close by the case we 
are considering for sufficiently small perturbations. We can 
construct a reference frame (r, 0, <f>) tilted midway between the 
two axes [9] , with the rationale that each semiperiod has its own 
adiabatic invariant, Eq. (6), tilted with respect to the other, and 
we expect that the composed trajectory will on average lie on 
the adiabatic invariant midway between the two. We can show 
that this will be the case in the limit when the axis separation 
and the Reynolds number are sufficiently small. The change 
in r and over a period is bounded by an quantity of 0(a), 
while 4> changes by an arbitrary amount over the same period. 
Thus for small axis separation a and low Reynolds numbers, r 
and are action variables, and <f> is an angle. Having this set of 
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approximate action-angle variables, we are able to calculate 
the invariant surfaces and the resonant surfaces associated with 
a nonintegrable perturbation of a two-action flow. 



adiabatic invariant and the resonances for a particular value of 
the parameters. 
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FIG. 2. Level curves of the adiabatic invariant surfaces (thin lines) 
are shown inside the two spheres (thick lines) together with the low- 
est-order resonant surfaces intersecting them (thin dashed lines) on 
the plane y = r sin 8 sin</> = for a particular set of parameter values 
for small axis separation. The resonances shown are, from the inner 
sphere to the outer, k/n = -2/1, -3/2,-1/1,-1/2, 0/1, and 1/2. 

Although we are not able to obtain an exact expression for 
the adiabatic invariant of the stroboscopic map for all a , we 
can obtain perturbatively in a an expression 

Re f A n % < 

H = — + A 21 + A 3i r 3 + A 4l r 5 
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(7) 



valid for small axis separations, where a, ;j , and A, ;j represent 
the constants a; and A{ evaluated for the two semiperiods 
j = 1,2. We can similarly obtain the resonant surfaces: in the 
tilted coordinate system (r, 6, <f>), <f> correspondingly strobed 
is the fast angular variable. The resonant rotation of <f>, as in 
Liouvillian maps, occurs when = lirk/n, where k and n 
are integers. For small axis separations, this is true when 
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where T is the axis-switching period. Notice that these reso- 
nant surfaces are independent of 6, i.e., are spherical shells of 
varying radii r^/ n , up to first order in the axis separation. In 
Fig. 2 we show the locations in the meridian plane of both the 
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FIG. 3. (a) The projection of a strobed trajectory onto the plane 
y = r sin 8 sin <f> = 0. The parameter values were chosen to be the 
same as in Fig. 2, such that the lowest-order 0/1 resonance (dashed) 
is in the middle of the region between the spheres, while leaving other 
primary resonances distant. The axis separation angle is a = 0.1° 
and the Reynolds number Re = 0.1. (b) Time evolution of the 
adiabatic invariant H for this case. 

Numerically computed strobed trajectories show that the 
theoretical picture we painted above is accurate and that the 
phenomenon of resonance-induced diffusion, previously re- 
ported only for Liouvillian maps, is present here with strik- 
ing resemblance in a real fluid flow. To illustrate this, we 
choose the flow parameters such that the relative geometry 
of the adiabatic invariants and the resonances matches that of 
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Fig. 1(a) above: the lowest-order resonance in Fig. 2 passes 
near the middle of the family of adiabatic invariants. In 
Fig. 3(a) we show a projection on the meridian plane of the 
stroboscopically-sampled trajectory of the flow at these pa- 
rameter values. The similarity between Fig. 3(a) and Fig. 1(a) 
is immediately obvious. Furthermore, in Fig. 3(b), we show 
the strobed time evolution of the invariant given by Eq. (7). 
Again, its likeness to Fig. 1(b) is undeniable, which highlights 
both the presence of resonance-induced diffusion here, and the 
accuracy of our assumptions about the adiabatic invariant. 
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FIG. 4. A slice of a strobed trajectory between— 0.01 < y < 0.01 
for the same parameter values as Fig. 3. Locations of low-order reso- 
nances are shown as in Fig. 2. Forty thousand points (corresponding 
to many times this number of periods) are plotted here, all from the 
same initial condition. 

Finally, as the slice of a strobed trajectory in Fig. 4 attests, 
the mechanism of resonance-induced diffusion is responsi- 
ble for the existence of asymptotically-globally-space-filling 
trajectories. This slice illustrates how effective resonance- 
induced diffusion is at mixing; the perturbation, which is of 
0(a), is very small here, yet the only regions where the tra- 
jectory has not yet ventured are close to the poles, and the 
middle of the vortices, which lie inside the 0/1 resonance 
shell. In time, the trajectory will diffuse into these regions 
too, through the action of the higher-order resonances that are 
dense throughout the space. Increasing the axis separation a 
increases the amount of time the trajectory is captured into 
resonance, and with it the size of the jumps and the diffusion 
rate. Other quantities relevant to mixing — the stretching and 
folding measured by the Lyapunov exponents of the flow — 
also increase as powers of a. Experimental control of the 
mixing rate can then be achieved through adjustment of the 
strength and density of the resonances with the axis separation 
a, the angular velocity ratio of the spheres, and the period of 
the motion [9]. Further quantitative studies of mixing effi- 



ciency, as well as the effects of the interplay of space-filling 
trajectories and molecular diffusion on mixing and transport, 
are now in progress. 

In summary, we have shown that the properties of nearly- 
integrable two-action Liouvillian maps are highly relevant 
to the transport features of a large class of real fluid flows. 
We expect to see in the stroboscopic maps of such flows 
the resonance-induced diffusion characteristic of two-action 
maps, consisting of motion on invariant surfaces interspersed 
with periods of motion on resonant surfaces. We have also 
shown that our action-angle classification of Liouvillian maps 
has good predictive capability. This finding opens up an 
avenue for the experimental verification of the existence of 
space-filling trajectories in this and other similar flow geome- 
tries in the presence of small periodic modulations. There are 
many technological implications of this novel diffusion phe- 
nomenon, and we aim to have excited the reader about the new 
possibilities for enhanced mixing that may arise. 
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